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We investigate the effect of disordered vacancies on the normal-state electronic structure of the 
newly discovered alkali- intercalated iron selenide superconductors. To this end we use a recently de- 
veloped Wannier function based method to calculate from first principles the configuration-averaged 
spectral function (A(k,cj)) of K0.sFe1.6Se2 with disordered Fe and K vacancies. We find that the 
disorder can suppress the expected Fermi surface reconstruction without completely destroying the 
Fermi surface. More interestingly, the disorder effect raises the chemical potential significantly, giv- 
ing enlarged electron pockets almost identical to highly doped KFe2Se2, without adding carriers to 
the system. 



PACS numbers: 74.70.-b, 71.15.-m, 71.18.+y, 71.23.-k 

The newly discovered [1 alkali-intercalated iron se- 
lenide selenide superconductors A ;E Fe2- 2/ Se2, with A=K, 
Rb, Cs, share a set of properties that distinguish them 
from any of the other iron based superconductor (FeSC) 
families. Their transition temperature is found to be as 
high as T c = 31K at ambient pressure [1 j, high in com- 
parison with the other 11-type chalgogenides. They are 
the first FeSC on the verge of becoming an antiferromag- 
netic ( AFM) insulator [2] . The ( AFM) order has a novel 
block type structure with an unprecedented high transi- 
tion temperature of T/v = 559K and magnetic moment of 
3.31/i B /Fe [3 . From ARPES experiments [HHP], only 
electron pockets are found and no hole pockets, seemingly 
ruling out the popular s ± pairing scenario [TT]. Last, but 
from an electronic structure point of view certainly not 
least, they contain significant numbers (~ 20%) of Fe and 
K vacancies (Vp e and Vk)- 

While the presence of these vacancies has been well 
accepted, their distribution is currently under debate. 
From the superlattice peaks in neutron diffraction studies 
[3] it is concluded that the Fe vacancies form a y/E x y/E 
superlattice with strong long-range order. On the other 
hand, a different scenario has been suggested for those 
samples which display simultaneously long range mag- 
netic order and superconductivity, namely that of a phase 
separation, whereby the ordered Fe vacancies only ex- 
ist in one of the phases. Some scanning tunneling mi- 
croscopy [12] , nuclear magnetic resonance [13] , and trans- 
mission electron microscopy (TEM) measurements [14] 
conclude that the second phase contains no Fe vacancies. 
Other TEM [4j[15] and x-ray absorption near edge struc- 
ture measurements [16] conclude that the second phase 
contains disordered Fe vacancies. 

An obvious puzzle emerges immediately upon examin- 
ing angle-resolved photoemission spectroscopy (ARPES) 
measurements [4HT0]. As was demonstrated theoreti- 



cally [17 , the y/E x y/E superlattice of Fe vacancies must 
induce a strong reconstruction of the Fermi surface and 
the band structure, due to the huge scattering associ- 
ated with the vacancy. However, current ARPES results 
all show a simple electronic structure without any indi- 
cation of reconstruction. A possible solution [17] to this 
puzzle is that the regions to which ARPES measurements 
are sensitive must contain poorly ordered Fe vacancies. 
This explanation, however, is vulnerable to the potential 
problem that while vacancy disorder might suppress the 
band structure reconstruction, it may also destroy com- 
pletely the Fermi surface. It is thus timely to examine 
the physical effects of Fe vacancy disorder in this class of 
materials. 

The ARPES experiments also reveal a peculiar behav- 
ior of the A ;E Fe2- 2/ Se2 family. Regardless of the nominal 
value of x and the observed Fermi surfaces always con- 
tains large electron pockets in the zone corner, as if the 
system is strongly electron doped compared to the "245" 
parent compound A0.8Fe1.6Se2. In fact, the chemical po- 
tential is shifted so much that not only do the hole pock- 
ets at the zone center submerge entirely below the chem- 
ical potential, but also a small electron pocket appears 
instead. Consequently, these systems are commonly re- 
garded as strongly electron doped, closer to AiFe2Se2 
(122) than 245, without a reasonable account of the ori- 
gin of the large number of extra electrons. Considering 
the essential role of doping in the magnetism and super- 
conductivity of these materials, a physical explanation of 
this peculiar behavior is clearly of highest importance. 

In this Letter, we address these two important issues 
by studying the influence of vacancy disorder on the 
alkali-intercalated Fe selenides, using a recently devel- 
oped Wannier function based method to calculate from 
first principles the configuration-averaged spectral func- 
tion (A(k,(jj)) of K0.8Fe1.6Se2 with disordered Fe and K 
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vacancies. We find that the disorder does indeed suppress 
the band structure reconstruction without completely de- 
stroying the Fermi surface. Furthermore, we find that 
the disorder raises the chemical potential significantly, 
eliminating the hole pockets and enlarging the electron 
pockets, in excellent agreement with the ARPES exper- 
iments, without adding carriers. We provide a micro- 
scopic explanation for the origin of the effective doping 
in terms of the disorder induced frequency broadenings of 
the bands. In addition, we find the emergence of strongly 
incoherent carriers. Our result demonstrates clearly the 
general inapplicability of Luttinger theorem in disordered 
systems, and suggests an alternative interpretation of ex- 
isting ARPES measurements. 

The band structure of a disordered system is 
given by the configuration-averaged spectral function: 
(A n (k,u))=Y, c P c A c n (k,uj) of Wannier orbital [18] n, 
crystal momentum k and frequency cj, in which configu- 
ration c is weighted by its probability P c . By treating the 
disordered configurations within the supercell approxi- 
mation, their spectral functions A^(k, uo) can be obtained 
directly from the supercell eigenvectors and eigenvalues, 
using the unfolding method [19]. To handle the computa- 
tional cost related to the large sizes of the supercells, es- 
sential for a proper treatment of the disorder, we employ 
the recently developed Wannier function based effective 
Hamiltonian method [20]. The low energy Hilbert space 
is taken in [-6,2] eV consisting of the Wannier orbitals of 
Fe-d and Se-p characters. The influence of the Fe and 
K vacancies is extracted from three DFT [21] calcula- 
tions: the clean KFe2Se2 and the single vacancy super- 
cells K7Fei6Sei6 and KgFeisSeie- For the configurational 
average, we use 10 large supercells (e.g. Figjljb)) of ran- 
dom size, shape and orientation containing 176 atoms 
on average, and assume an equal probability among the 
configurations. Benchmarks demonstrating the high ac- 
curacy of the effective Hamiltonian method and the con- 
vergence of the configuration average against the size and 
the number of configurations are given in Ref. [21 . 

As a reference, let's first re-examine the fully y/E x y/E 
ordered 245 system that shows strong band structure 
reconstruction [17] . To facilitate the comparison with 
ARPES, Figure[TJc) shows the band structures and Fermi 
surfaces in the 1-Fe Brillouin zone (c.f. Fig. [TJa)). In- 
deed, the lower panel shows a very strong reconstruc- 
tion of the Fermi surface that bears no resemblance to 
the observed ARPES measurements [4lfT0 . This strong 
Fermi surface reconstruction is not a special property of 
the y/E x y/E ordering, but expected to happen for other 
vacancy orderings [I4j [15] in general, due to the very 
strong scattering against the Fe vacancies. The top panel 
also shows that the whole band structure, not just the 
Fermi surface, is reconstructed, with gap openings and 
complex folded features occurring everywhere in both Fe 
and Se bands. For example, the red bars in the middle 
panel indicate two large gap opening on the order of 300- 



400 meV. These kinds of strong reconstructions should in 
principle allow ARPES to quantify the strength of the va- 
cancy order, including the semiconducting samples that 
exhibit no Fermi surface. The fact that no such recon- 
struction has been reported is at odds with the strong 
vacancy order observed by neutron scattering [3]. 

The configuration- averaged spectral function in 
Fig. [ljd) confirms that strong disorder in the vacan- 
cies introduces very large incoherent scatterings. The 
Fe bands now develop large linewidths in both momen- 
tum and frequency, reflecting their short mean free path 
and lifetime. More dramatically, the top of the Se bands 
within [-4,-2] eV even becomes unrecognizable, due to the 
large repulsive potential of ~0.6eV/0.2eV from Fe/K va- 
cancies on their 4/8 neighboring Se orbitals. Neverthe- 
less, even for such a completely disordered case, Fig.JIJd) 
shows that the Fe band structure and the Fermi surface 
survive the strong vacancy disorder effects, displaying the 
characteristics of a dispersive band structure. On the 
other hand, Fig. [ljd) shows the clear suppression of the 
reconstruction in the band structure and the Fermi sur- 
face. All the gap openings and folded features are washed 
out. In fact, other than the strong broadening effects and 
a clear reduction of bandwidth, the resulting band struc- 
ture and Fermi surfaces agree very well with those of the 
clean 122 system in Fig. [lje). We have thus provided 
quantitative support to the notion that the regions to 
which ARPES measurements are sensitive must contain 
poorly ordered Fe vacancies. 

Figure [ljd) also shows an unexpected effect of the va- 
cancy disorder, namely a large shift of the chemical po- 
tential. Indeed, the chemical potential is now ~200meV 
above the tip of the hole pockets, and allows a small elec- 
tron pocket to appear around the Z point. Correspond- 
ingly, the electron pockets at the X/R point are much 
larger than the parent compounds of all families of Fe- 
based superconductors. In fact, the resulting Fermi sur- 
face looks exactly like that of the heavily doped 122 sys- 
tem in Fig.[TJe), in excellent agreement with the ARPES 
measurements, even though the system actually contains 
no additional electrons. (A Luttinger count of the car- 
rier density from the volume of the Fermi surface would 
give the illusion of -f 0.5/Fe doping for this undoped sys- 
tem.) This is a clear demonstration that Luttinger theo- 
rem should never be applied to systems with strong disor- 
der [22, 23 . Specifically in this case, our result suggests 
that the experimental samples showing 122 like Fermi 
surface might not even be doped at all. 

Microscopically, this strong shift of chemical potential 
(or an effective "doping") is to be distinguished from 
real doping, for example in the case of Co substituted 
BaFe2As2, in which the chemical potential shift is a di- 
rect consequence of the additional valence electron in the 
system [23]. (This conclusion differs from the previous 
studies [24] [25]). Figure [2|a) reveals more insight by ig- 
noring the impurity potential in the calculation (but still 
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FIG. 1: (color online) (a) An example of a large sized supercell used for the conngurational average and (b) 1-Fe Brillouin 
zone. Band structure, density of states and Fermi surface of v^x ordered Fe and K vacancies(c), disordered Ko.8Fei.6Se2(d) 
and KFe2Se2(d). The dashed lines in the Fermi surface plots mark the 2-Fe Brillouin zone boundaries. The band structure and 
Fermi surface of the disordered K0.sFe1.6Se2 are enhanced by a factor of 2 and 10 respectively for better visualization. 
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FIG. 2: (color online) Spectral function obtained from dis- 
ordered removal of Fe orbitals alone (without the impurity 
potentials) (a), and from the virtual crystal of the impurity 
potentials (without the removal of Fe orbitals). Both show 
negligible effective doping in comparison with the full disor- 
dered calculation. (c)(d) Illustration of microscopic mecha- 
nisms contributing to effective doping through broadening in 
frequency. 



removing Fe orbitals at the disordered vacancy sites). 
In this case, the effective doping becomes negligible, as 
an "expected" chemical potential of an undoped system 



is found, giving both hole and electron pockets of reg- 
ular sizes. On the other hand, an average account of 
the impurity potential via the virtual crystal approxima- 
tion (c.f. Fig. [2^b)) misses entirely the effective doping 
as well. Therefore, the strong chemical potential shift is 
only realized from the disordered nature of the impurity 
potential. 

Figure [2] illustrates two microscopic mechanisms that 
contribute to the effective doping through disorder in- 
duced broadenings in frequency. For simplicity, let's con- 
sider the broadening of the momentum independent spec- 
tral function A(uj). First, in Fig. [2|c), when the spectral 
weight near the chemical potential (within the frequency 
range of the broadening) is larger below, a symmetric 
broadening in frequency pushes the chemical potential 
up. Second, in Fig. |2|d), an asymmetric broadening 
skewed toward higher frequency alone does the same[26j. 
The above test cases in Fig. [2|a) and (b) suggest that 
the essential mechanism in K ;E Fe2- 2/ Se2 is the aymmetric 
broadening resulting from the repulsive impurity poten- 
tial. 

Without increasing the total electronic density of the 
system, the significant shift of chemical potential leads 
to many important physical consequences through, for 
example, the change of Fermi velocity and the Fermi 
wave vector, in a manner similar to real doping in regular 
systems without disorder. In particular, the significant 
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FIG. 3: (color online) Disorder induced spectral weight 
distributions, represented in the 2-Fe Brillouin zone, (a) as 
a function of crystal momentum k and orbital indices n, 
and at fixed crystal momenta ko and band index j with (b) 
ko = 0.8ZR and j = 12 belonging to the large electron pocket 
in the zone corner, (c) ko = 0.5ZR and j = 10 belonging to 
the hole pocket leaking into the Fermi surface, (d) ko = 0.5ZR 
and j = 14 displaying an incoherent satellite at the Fermi sur- 
face and (e) their integrations displaying ^20% loss of spectral 
weight due to the Fe vacancies. 



suppression of phase space for long-range magnetic order 
around (7r,0) resembles very much the effect of real dop- 
ing. This offers a natural explanation of the occurrence 
of a superconducting phase in this undoped system, and 
suggests a novel mechanism of promoting superconduc- 
tivity via strongly disordered impurity potentials. 

Interestingly, there have indeed been experimental ob- 
servations suggesting enhancement of superconductiv- 
ity via disorder in these systems. Reference [27] re- 
ported that the same sample can be made superconduct- 
ing/magnetically ordered through quenching/ annealing. 
High pressure experiments [28] reported a seemingly sec- 
ond superconducting transition at considerably higher 
temperature T c = 4SK at pressure above 11.5GPa, con- 
sistent with promotion of superconductivity via pressure- 
induced suppression of vacancy order. A second super- 
conducting phase has also been recently reported at am- 
bient pressure [29] . 

Obviously, for disordered impurities to promote super- 
conductivity, their positive influence (from the above ef- 
fective doping or other disorder effects discussed below) 
must overcome the negative influence of disorder (for 
example pair-breaking or phase fluctuation scattering). 
This issue is intimately tied to the nature of the scatterer 
(for example, magnetic or nonmagnetic, weak scattering 
or strong scattering regime, and so on) and the super- 
conductivity itself (for example, sign changing order pa- 
rameter or not, dominant amplitude fluctuation or phase 
fluctuation, and so on). We anticipate rich physics to 
be revealed from future quantitative investigations along 
this direction. 

It is also important to account for the removal of Fe 
orbitals in the vacancy sites. Indeed, representing the 
spectral function in the eigenbasis \kj) of crystal mo- 
mentum k and band index j of the uniform 122 system, 



Aj(k,uS) = -±lm{kj\G(u)\kj) [23 , Fig. [5] shows that the 
total spectral weights of all well-defined low-energy quasi- 
particle drop to ~0.8, in good accordance to the 20% 
Fe vacancy. As advocated previously [23] , such a large 
reduction in the spectral weight is expected to change 
strongly (even qualitatively) the standard picture of spin 
fluctuations [30] [31], through a significant reduction of 
the bare susceptibility by roughly a factor of 0.8 2 ~ 0.6. 
The general argument [23] of suppression/resilience of 
phase space of antiferromagnetic/superconducting order 
should remain valid, which leads to a potential promotion 
of superconductivity. 

In addition, our result reveals strong incoherent fea- 
tures in this system. For example, Fig. |3|c) illustrates 
that small amount of diffusive hole carriers in the sys- 
tem are still present in this seemingly highly electron 
doped system, through the long tails of the hole bands. 
This suggest that hole carriers could still contribute to 
the effective pairing [31]. Similarly, Fig. [3]^d) shows the 
presence of highly incoherent carriers right at the chem- 
ical potential, originating from the normally unoccupied 
band 0.5eV above the chemical potential. All of these 
strong effects illustrate vividly the rich physics induced 
by disorder, well beyond a simple modification of the 
band filling or dispersion. 

Finally, it is important to clarify that our case study 
with fully disordered vacancies cannot correspond exactly 
to the real materials, since that would not be consistent 
with the observed superstructures from diffraction exper- 
iments [3j HJ [I4j [15] . It is instead meant to illustrate the 
general features of poorly ordered regions of the sample, 
for example, near the surface where diffusion of vacancies 
is easier, or in the vacancy disordered region of a phase 
separated sample. Also, the vacancies are likely to be 
short range ordered (SRO), given their charged nature. 
Furthermore, additional scattering related to short-range 
magnetic and orbital correlation [32j [33] should also be 
present. While we do not expect such SROs to quali- 
tatively change our conclusions, further studies on their 
effects are desirable. 

In conclusion, we studied from first-principles the 
normal state configuration-averaged spectral function 
(A(k,uj)} of K0.8Fe1.6Se2 with disordered vacancies. We 
found that strong vacancy disorder does not destroy the 
band structure and Fermi surface, but is able to suppress 
the significant band structure reconstruction expected 
in vacancy ordered systems, thus producing the smooth 
quasi-particle band structure observed by ARPES. Fur- 
thermore, we found an intriguing disorder-induced "ef- 
fective doping" that enlarges significantly the electron 
pockets to be comparable to the highly electron doped 
KFe2Se2, in excellent agreement with ARPES observa- 
tions. Our result demonstrates clearly the general inap- 
plicability of Luttinger theorem in disordered systems, 
and suggests that the region to which ARPES measure- 
ments are sensitive might be almost undoped after all, 
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but with strong vacancy disorder. Finally, the observed 
effective doping, together with a strong (~20%) loss of 
coherent spectral weight, offers a novel, general mech- 
anism to promote superconductivity in the presence of 
strong impurity potentials, as indicated by several cur- 
rent experiments. 

Work funded by the U S Department of Energy, Office 
of Basic Energy Sciences DE-AC02-98CH10886 under the 
CMCSN program. 
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DETAILS OF THE DENSITY FUNCTIONAL THEORY CALCULATIONS 
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FIG. 4: (Left) normal cell KFe2Se2, (middle) single Fe vacancy supercell KsFeisSeie and (right) single K vacancy supercell 
K 7 Fei 5 Sei 6 . 

The lattice structure of the normal cell, KFe2Se2, is depicted in the left of figure [4] Its space group I4/mmm and its 
lattice parameters a=7.46 bohr and c=26.76 bohr and z=0.354, were taken from [2]7 The primitive normal cell lattice 
vectors expressed in Cartesian coordinates are given by: a\ — — ^ax+^ay+^cz ; CL2 = \ax — \ay + \cz ; as = 
\ax + \ay — \cz . To capture the influence of the Fe and K vacancies an eight times larger supercell is used to properly 
treat the non-local influence on the nearest As and Fe orbitals. The corresponding KgFeisAsie/K^FeieSeie supercells 
are depicted on the right / middle of figure [4] and its primitive lattice vectors expressed in normal cell lattice vectors 
are given by: A\ = 2<22 + 2*23 ; A2 = 2a\ + 2*23 ; As = a\ + a<i. We applied the WIEN2K[1 implementation of the 
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full potential linearized augmented plane wave method in the local density approximation. The k-point mesh was 
taken to be 10x10x10 for the undoped normal cell and 5x5x3 for the supercells respectively. The basis set sizes 
were determined by RKmax=7. 



BENCHMARKS OF THE EFFECTIVE HAMILTONIAN AGAINST DFT 



To explore the accuracy and efficiency of the effective Hamiltonian method [3] for the case of K :c Fe2- y Se2, we 
present comparisons of spectral functions A n (k,uj) calculated from the full DFT and the effective Hamiltonian (see 
figures [5][7]) . The size of the deviations between the full DFT and the effective Hamiltonian should be compared with 
the size of the impurity induced changes. For this purpose the spectral function of the undoped KFe2Se2 is also 
plotted as a reference for each benchmark. The test cases in figures [5] and [6] are designed as extreme test cases for 
the "linearity" approximation of the impurity influence (see formula (1) of Ref. [3 ). For testing the linearity of the K 
vacancy influences, the case of Fe2Se2 corresponds to the maximum extrapolation from (1 — x) = 0.125 in the single 
impurity cell (see right of Fig. [4]) to (1 — x) = 1. For testing the linearity of the Fe vacancy influences, a smaller 
extrapolation is more meaningful, in order to have remaining Fe bands. The bands of KFeSe2 are very reasonably 
described, despite the extreme extrapolation in doping from y = 0.125 in the single impurity cell (see middle of Fig. 
[I| to y = 1. The test case in figure [7] of a y/5 x y/5 supercell is designed to test the partitioning of the impurity 
influence from its super images (see section IV of the supplementary information of Ref. [3]). The basis set of Linear 
Augmented Plane Waves (LAPW's) used in the full DFT is ~ 30 times larger then the basis set of Wannier functions 
used in the effective Hamiltonian method. Since the number of floating point operations of diagonalization depends 
cubically on the size of the matrix this implies an efficiency increase by a factor of 30 3 « 3 • 10 4 . Furthermore the 
full DFT calculation involves multiple self consistent cycles (typically 15) whereas the effective Hamiltonian method 
(as currently implemented) requires only a single diagonalization, which increases the efficiency by another order of 
magnitude to ~ 6 • 10 5 . 
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CONVERGENCE WITH RESPECT TO SIZE AND NUMBER OF CONFIGURATIONS 



10 configurations of 176 atoms on average 20 configurations of 352 atoms on average 
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FIG. 8: 

In figure [8] we demonstrate the convergence of the spectral function of disordered K0.sFe1.6Se2 with respect to the 
number and the size of the configurations. 
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